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Abstract: The recently proposed construction of chiral fermions on lattices with bound- 
aries is tested in an interacting theory up to first order of perturbation theory. We confirm 
that, in the bulk of the lattice, the chiral Ward identities take their continuum value up to 
cutoff effects without any tuning. Universal quantities are defined that have an expansion 
in the renormalised couplings with coefficients that are functions of the physical size and 
the periodicity in the spatial direction. These coefficient functions have to be identical 
for different discretisations. We find agreement with the standard Wilson fermions. The 
computation is done in the asymptotically free Gross-Neveu model with continuous chiral 
symmetry. 
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1. Introduction 

We present a 1-loop perturbative calculation that confirms the favourable properties of 
chiral fermions on a lattice with boundaries [|J . The calculation is performed in the chiral 
Gross-Neveu model Q. Before I go into the details and results of the calculation let me 
briefly introduce field theories on lattices with boundaries. 

In the Schrodinger functional (SF) of a quantum field theory the fields are defined on 
a d + 1 dimensional cylinder. The fields are subject to Dirichlet boundary conditions in 
the time direction (P± = |(1 ± 70)) 

P+^(x) = 0, ^(x)P_=0 at x = 0, (1.1) 

P-ip(x) = 0, lp(x)P + = at x = T (1.2) 
and periodic boundary conditions in the space directions 

ifj(x + Lk) = e ie ip(x) , 6»€[0,2vr), k = l,...d. (1.3) 

Discretising such a theory on a space-time lattice with boundaries has advantages. If 
the temporal extension is a multiple of the spatial extension, say T = 2L, the inverse of 
the spatial extension 1/L provides an infrared cutoff and is the natural energy scale in 
the massless theory. In fact this has enabled a fully non-perturbative determination of the 
scale dependence of the fundamental parameters of iVf = 2 QCD, see Q and references 
therein. 

However, naively one would expect that theories become more complicated on lat- 
tices with boundaries. First of all the distinction between different classes of boundary 
conditions, for example Dirichlet and Neuberger boundary conditions, really only makes 
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sense for smooth fields, which we do not have on the lattice. The boundary conditions 
therefore must be encoded in the lattice action and will arise dynamically in the contin- 
uum limit. In general this might involve the necessity to fine-tune some parameter in the 
action. Secondly the boundaries might cause additional divergences and thus lead to ad- 
ditional counter-terms. Both of these issues have been addressed for QCD in Q with the 
result that there is no need for fine-tuning and that one just has to add a renormalisation 
factor for the boundary fields. 

The fact that the SF boundary conditions arise naturally (without fine-tuning) in the 
continuum limit has been understood on the basis of dimensional counting and boundary 
critical phenomena, see [Q and references therein. The SF boundary conditions are part 
of the definition of the continuum limit. Together with the symmetry properties and the 
dimensionality of the system they form a SF universality class: any local discretisation has 
this continuum limit. 1 

Exactly this observation led the author of [P to a formulation of the lattice SF of QCD 
with chiral fermions. In the continuum the SF boundary conditions break chiral symmetry. 
Therefore lattice fermions that fulfil the Ginsparg- Wilson relation || on the whole lattice 
(with boundaries) can not have the right continuum limit. The minimal modification to the 
Ginsparg- Wilson relation that breaks the lattice chiral symmetry at the boundary would 
be a term A# that is supported at the boundary with exponentially decaying tails 

75D + D75 = a D~f 5 D + A B . (1.4) 

Given a local solution to this equation chiral Ward identities are expected to take their 
continuum form far away from the boundaries. 

As mentioned above, in principle there are many possible local lattice formulations 
that have the same continuum limit. In Q the author gives one solution to the modified 
Ginsparg- Wilson relation, a modified Neuberger-Dirac operator 

£, = s{ 1 - 1(^ + 75^75)} , (1-5) 

U = A (A^A + caP)" 1 ' 2 , a = — ^— , A = l + s — aD w , (1.6) 

1 + s 

where -D\v is the massless Wilson-Dirac operator in the presence of the boundaries as 
introduced in i.e. the standard infinite lattice Wilson-Dirac operator in the range 
< xq < T and at all other times the target field x = D\yip is set to zero. 2 The projector 
P is defined through 

P^x) = J{^P-^)|^ + W+^)U. B } , (1-7) 
and the parameters c and s must be chosen such as to ensure locality (see next section). 



The operator (|1.5|) solves (1.4) if the lattice spacing a is replaced by the rescaled value a. 



For all unexplained notation we refer the reader to (7j. 

1 Note that in [H] a SF is proposed with chirally rotated boundary conditions, which however break parity 



and therefore are distinct from the boundary conditions (1.1) considered here. 

2 This Dirac operator maps the space of fields defined at all xo, but set to zero at xo < a and xq > T - 
into itself. 
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2. Free fermions 



For Wilson fermions the free propagator in presence of the boundaries can be calculated 
explicitly and is given in ||] . In the case of chiral fermions we do not even have an analytic 
expression for the Dirac operator. To obtain the free propagator in closed form like in the 
Wilson case might be possible but is much more difficult. Nevertheless, the operator under 
the square root in ( |1.5[) can be worked out in the time-momentum representation. 

Taking into account the definition of _D\y i n the presence of the boundaries given above, 
the operator under the square root explicitly reads 

A*A + caP = (1 + s) 2 + sa 2 ^ d^d„ + ±a 4 ^ d%d„dtd v + (c - l)aP , (2.1) 

with forward d^ix) = ^(if)(x + afi) — ip(x)) and backward d*^{x) = ^(ip(x) — i/)(x — aft)) 
finite differences. It is Hermitian and therefore has real eigenvalues. Due to translation 
invariance in space the spatial eigenfunctions are plane waves with momentum values that 
are integer multiples oi2n/L in the range —ir/a < p^ < 7r/a. 3 For c = 1 the operator ( |2.1[ ) is 
diagonal in Dirac space, so its eigenvalues are 2 ( - d+1 ) //2 -fold degenerated. Its eigenfunctions 
(in space-time) are then given by 

T17T 

c = l: v POjP (x) = e tpx sm(p x ) , p = — , n = 1, 2, . . . ,T/a - 1 . (2.2) 



For c / 1 the operator fl2,l|) is not diagonal in Dirac space. The eigenfunctions are still 
given in terms of sin-functions 

c ^ 1 : w po , p (x) = e ipx {P_ sin (p x + b(p )) + P+ sin (p (T - x ) + b(p )) } , (2.3) 

where 

6(p ) = -arctan^3-^^l_j , q = ^p 2 - s , p, = § sin(a^/2) . (2.4) 

But the allowed values of po are now given by the solutions of the transcendental equation 

b(p ) = -p T. (2.5) 
In either case the eigenvalues are given by 

A POlP = (1 + s) 2 - sa 2 f + £ ^p 2 p 2 ■ ( 2 ' 6 ) 



Rewriting this as 



A P0 ,p = q a 2 p 2 + (1 + s) 2 - sa 2 p 2 + ^ Yl ' ^ 

k<l 



3 Greek indices fj,,v, . . . run from to d and Latin indices k, I, . . . from 1 to d. 
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one can easily show that the eigenvalues are bounded from below by (1 — |s|) 2 for \s\ < 1 
and c > 1. In particular, the constraint on c ensures that the combination qa 2 p^ is always 
positive. 

Adopting the argument in ref. |j, using expansion in Legendre polynomials, we con- 
clude that in the free theory the locality of the Dirac operator (|ll]) is guaranteed for this 
range of parameter values. 

The eigenfunctions ( [2.2; ) (or Q2.3| )) may be orthonormalised and used to write down an 
analytical expression for the kernel D(x,y) of the Dirac operator. But the evaluation of 
D(x, y) in this way would be very expensive, since it involves a sum over momenta po which 
in turn are determined for each set of parameter values by the roots of a transcendental 
equation. 

In our computation we used a different approach. First the operator ( |1.5|) is half- 
Fourier-transformed to its time-momentum representation D(xq,p). The square-root of 
the remaining (d+ 1)(T — 1) x (d+ 1)(T— l)-matrices (one for each value of p) is computed 
with the min-max polynomial as explained in section 4 of [10]. For the final step, to obtain 
the propagator, we used the built-in inversion routine of MATLAB. 

The so calculated modified Neuberger-Dirac operator has been inserted in the Ginsparg- 
Wilson relation to compute the difference Ab and to confirm, that it is in fact localised at 
the boundaries with exponentially decaying tails. 



3. Chiral Gross-Neveu model 

To test the operator Ql.5| ) beyond free fermions we introduce a two-dimensional theory 
with four-fermion interactions. The euclidean continuum action of the chiral Gross-Neveu 
model can be given in the form 

S8gn = / dx 2 " b 2 (O SS - Opp) - \g 2 v O vv } , (3.1) 

O SS = (^) 2 , Opp = (?75^) 2 , O vv = ^(V^Vf • (3.2) 

All Dirac and flavour indices are suppressed and contracted in a straightforward way. 
For N flavours of fermions this action possesses an [/(A)-flavour symmetry and an U{\) 
continuous chiral symmetry. For a detailed derivation of the possible terms (in terms of 
renormalisability) see chapter 5 of ref. Q . The terms in (|3.1| ) are a full set of allowed terms 
that respect the above mentioned symmetries (plus euclidean symmetry). For N > 2 the 
model shares with QCD the property of an asymptotically free coupling (namely g 2 ). 

Discretisation with Wilson fermions is straightforward. But since the Wilson term 
breaks chiral symmetry a mass term and an additional coupling have to be added 

Scgn.w = ffl2 ^2 ( Dw + m °) ^ ~ \9w{Oss - Opp) - \S 2 P Opp - \gv. w Ovv) ■ (3.3) 

X 

No additional coupling (or mass term) is necessary if one uses a Dirac operator that 
solves the Ginsparg- Wilson relation. Since such an operator comes with a lattice chiral 
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symmetry [11|, there are no other allowed dimension 2 operators. However, to make the 
action manifestly invariant under that symmetry, one has to add irrelevant dimension 3 
operators via the substitution tp —> t[? = (1 — ^D)tp in the four-fermion interaction terms. 
The action then reads 

SCGN.GW = « 2 D ^ ~ hdgw(0ss ~ Opp) ~ ^9v,gwOvV^ ■ (3.4) 

X 

Defining now the SF of the Gross-Neveu model is straightforward. Since it is asymptot- 
ically free (for N > 2, which we assume from now on), the scaling dimension of local fields 
is equal to their engineering dimension. Therefore the argument about the naturalness of 
the SF boundary conditions in QCD Q holds also in the chiral Gross-Neveu model. And 
in [|7|] it is shown that just one additional renormalisation factor for the boundary fields 
has to be added. 

4. Chiral Ward identity 

The continuum chiral Ward identity 

(d lx A fl (x)0(y))=2m(P(x)0(y)), x^y (4.1) 

implies on the lattice for vanishing renormalised mass m,R = 

{{0) K d,{A il ) K (x)) = 0(a). (4.2) 

This can be used to compute the critical value of the bare mass parameter mo and the chiral 
symmetry restoring value of the bare coupling 5p in fl3,3| ) to second order in perturbation 
theory (PT) §. 

As indicated in (|4.2| ) there are chiral symmetry breaking effects that are lattice arti- 
facts. In this section we compute these effects for Wilson fermions and the proposed chiral 



fermions. To this end correlation functions that match the operators in (4T) and (^ 
have to be defined. 

Correlation functions like 

fA(x ) = ~ Yl (#07075^) C(Vi)T5C(*0 > (4-3) 

2/1,21 

correlate boundary states (here a pseudo-scalar) with current insertions (here axial current). 
The boundary fields ((x±) = P-ip(x) and ((xi) = ip(x)P+ at xq = a become the non- 
vanishing components of ip, i/j at boundary xq = in the continuum limit. Similarly 
fp(xo) is defined to correlate a pseudo-scalar boundary state to a pseudo-scalar density 
insertion. (All correlation functions will be evaluated for Wilson and for chiral fermions. 
In the later case the fermion fields, bulk and boundary, are substituted as indicated above 
eq. 0.) 

With the help of these two correlation functions the bare current mass can be computed 
in PT 

r(xd) = t°/f (X f = + r ^)9 2 + 0(5 4 ) (4.4) 
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Figure 1: Bare current mass at 1-loop in PT (left). Logarithmic plot of this quantity for chiral 
fermions (right). The plots are for a 16 x 32 lattice, 9 = 0, N = 2 and g% = 0. 



If the renormalised mass is set to zero {tur = 0) by demanding Ldo (/a(^o))r = 0(a/L) 
for all xq and 9, this quantity is a direct measure of the lattice artifacts on the right hand 
side of (gj). 

At tree-level ro(xo) is 0((a/L) 2 ) for all xq for Wilson and for chiral fermions. At 1-loop 
the behaviour is similar for both discretisaions. There are chiral symmetry breaking effects 
localised at the boundaries, that survive the continuum limit and decay exponentially with 
the distance to the boundaries, see fig. [j]. 



5. Renormalised coupling 

Having confirmed the desired chiral properties of the proposed chiral fermions in the SF, 
we now renormalise the couplings at vanishing renormalised mass and define universal 
quantities that have a finite and unique continuum limit. 

To this end we define boundary to boundary correlation functions 

A = - 2 ( a/-V £ <C , («i)75A a C / (^i)C(yi)75A a C(^i)> (5.1) 

uiviyizi 

/2 = -^E<CK)C / (^i)), (/ 2 ) R = ^ 2 /2 and (/ 4 ) R = Z c 4 / 4 (5.2) 

U\Z\ 

and the renormalisation factor free ratio 

W = 77^- 1 = 77^- 1 ^ m R = (5.3) 

which depends on the periodicity angle 9 in the spatial directions. The A a s are the gener- 
ators of the algebra of SU(N) and C, £ are the equivalents of (, ( at xq = T — a. 

Renormalised couplings can now be defined (up to normalisation factors) as the dif- 
ference of R(9) at two values of 9 

f oc R[9{) - R(9 2 ) (5.4) 
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and similar for g v @. (The normalisation factors are functions of 0\ and #2 only and 
chosen such that in the continuum limit g 2 ~ g 2 for small g 2 .) For definiteness we fix one 
of them to 6\ = 0. There are then renormahsed couplings g 2 , c/y^ for each value 62 = f]. 
They have an expansion in PT and the 1-loop result for g 2 with the proper normalisation 
is 

f v °=° dl + gi b Ha/L) + c w g { V ) + 0(<7 6 ) (5.5) 
for Wilson fermions and 

-2 a^O 



g- v = 9 2 gw +gUboHa/L) + c g g w (rj) + 0{g 6 ) 



(5.6) 



for chiral (gw) fermions. In this formulae, 60 = —N/ir, is the correct universal first co- 
efficient of the beta- function (showing asymptotic freedom), which was computed earlier 
in the continuum The finite parts found in the two computations differ as expected. 
The result for the renormahsed couplings g v , corresponding to the vector- vector interac- 
tion, can be found in ref. |?]]. We here only note, that after rearranging the terms in the 
action with the help of Fierz-transformations, this coupling has (up to 1-loop) a vanishing 
beta-function. 

An universal renormalisation group invariant (RGI) quantity can now be defined as 
the difference of the renormahsed coupling for two different values of 77. This is because in 



the continuum a non-zero 9 in (1.3) shifts the momenta of the external legs, but does not 
effect the loop integrals. The difference 



FnGi(ji,a/ L ) = g v ( a / L ) - g V0 ( a / L ) > % = 1 • 

has an expansion in the renormahsed couplings (we omit the subscript rf) 

F KG1 (r,,a/L) = F?\r,,a/L)? + ( V ,a/L) g 2 g 2 v + (r,,a/L) g v + 0(g e ) 



(5.7) 



(5.8) 



where the coefficients are given by the finite parts. Since it is an universal quantity, these 
coefficients have a finite and unique continuum limit independent of the discretisation. The 
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Table 1: First order coefficients in the expansion of Frgi with Wilson (w) and chiral (gw) fermions. 

continuum values for different choices of 7/ are given in tab. |. The errors are estimated 
by fitting the first few terms of the Symanzik expansion of lattice diagrams to the data 
at lattice sizes from 4 x 8 to 64 x 128 (128 x 256) for chiral (Wilson) fermions (method 
explained in appendix D of |l3|]). The continuum values agree within the error, thus 
confirming universality for the operator (|1.5|). 



6. Final remarks 

The formulation of the SF for fermionic models of the Gross-Neveu type in has en- 
abled the first test of the recently proposed chiral Dirac operator in the presence of the 
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SF boundary conditions. At 1-loop we have shown that a chiral Ward identity takes its 
continuum value far away from the boundaries. The calculation of an universal quantity 
gives the right discretisation independent value. 

The size of the lattice artifacts has not been discussed so far. What we observe is that 
the chiral fermions are tree- level 0(a)-improved, if the value of c is tuned correctly. At 
1-loop it is not enough to tune c, because there are dimension 3 operators at the boundary, 
that spoil automatic 0(a)-improvement. This however is peculiar to the studied kind of 
models with four-fermion interactions. 

As next step, it would be desirable to perform a similar (perturbative) computation 
in QCD in order to provide guidance for the final non-perturbative application of chiral 
fermions in the Schrodinger functional. 
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